Trends in graphene edge properties and flake shapes: a first-principles study 
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We use density functional theory to determine the equilibrium shape of graphene flakes, through 
the calculation of the edge orientation dependence of the edge energy and edge stress of graphene 
nanoribbons. The edge energy is a nearly linear function of edge orientation angle; increasing 
from the armchair orientation to the zigzag orientation. Reconstruction of the zigzag edge lowers its 
energy to less than that of the armchair edge. The edge stress for all edge orientations is compressive, 
however, reconstruction of the zigzag edge reduces this edge stress to near zero. Hydrogen adsorption 
is favorable for all edge orientations; dramatically lowering all edge energies and all edge stresses. It 
also removes the reconstruction of the zigzag edge. Using the new edge energy data, we determine the 
equilibrium shape of a graphene sheet (with unreconstructed edges) to be hexagonal with straight 
armchair edges in the presence and absence of hydrogen. However, zigzag edge reconstruction 
produces graphene flakes with a six- fold symmetry, but with rounded edges. This shape is dominated 
by near zigzag edges. The compressive edge stresses will lead to edge buckling (out-of-the-plane of 
the graphene sheet) for all edge orientations, in the absence of hydrogen. Exposing the graphene 
flake to hydrogen dramatically decreases the buckling amplitude. 
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I. INTRODUCTION 

Graphene, a one-atom thick single sheet of sp 2 - 
bonded carbon atoms in a honeycomb lattice, is a zero- 
gap semiconductor (semi-metal) that exhibits extraor- 
dinarily high electron mobility and shows considerable 
promise for applications in electronic and optical de- 
vices, high sensitivity gas detection, ultracapacitors and 
biodevices^'^^^^ ' 10 : 11 . Graphene can be separated 
from bulk graphene using a peeling technique known as 
mechanical exfoliation^. Lithography techniques^ 2 - have 
recently been applied to fabricate patterned graphene 
nanoribbons (GNR) with widths below 100 A; such GNR 
were used to produce graphene field effect transistors 
with on-off ratios as high as ~ 10 7 . There has been 
concerted efforts to investigate the stability of graphene 
sheet o 13 ' 14 ! 15 ! 16 : 17 ! 18 . Since applications of graphene nec- 
essarily employ sheets of finite extent, it is also of interest 
to inquire as to the equilibrium shape of such finite sheets 
or flakes. 

Two distinct definitions of shape arise in the context 
of graphene sheets. The first is associated with a finite, 
perfectly flat graphene sheet. For such a sheet, the equi- 
librium shape is determined primarily by the energy of 
the graphene sheet edges. Edge energies are related to 
the thermodynamically stable shape in exactly the same 
way that surface energies are used to predict the equi- 
librium shape of a three dimensional crystal through the 
Wulff construction^. The second is associated with equi- 
librium elastic distortions of the graphene sheet as a re- 
sult of edge stresses (analogous to surface stresses for 
a three-dimensional crystal). Several recent studie a 15 : 16 
have shown that the two highest symmetry graphene 
edges have compressive edge stresses, leading to clastic 
distortions of the regions near these sheet edges (edge 



warping). In this paper, we employ density- functional 
theory (DFT) to determine graphene edge energies and 
edge stresses as a function of edge orientation both in 
the absence and in the presence of hydrogen atoms. We 
employ these results to investigate the equilibrium shape 
of a finite graphene sheet; both for the flat sheet (a two- 
dimensional intrinsic manifold) and one for which elastic 
distortions are permitted (a two-dimensional sheet em- 
bedded in three dimensions). 

II. CALCULATION METHODS 

In order to determine the properties of the edges 
of graphene sheets, we employ density-functional the- 
ory (DFT) 20 ' 21 within the local density approximation 
(LDA) using a self-consistent pseudopotential method, as 
implemented in the Siesta code^ 2 -. A double-^ plus po- 
larization basis set is used for the localized basis orbitals 
along with a large energy cutoff (400 Ry) for the real 
space integration. An electronic temperature of 0.01 eV 
is used to smear the Fermi-Dirac distribution. In order 
to study graphene sheet edges, we focus on flat graphene 
nanoribbons (GNR) that are periodic in one direction 
and are separated by 16 A of vacuum in the other two 
directions. In the direction parallel to the ribbon axis, 
we use a fc-point sampling criterion such that the number 
of k points rik in the periodic direction is determined by 
the smallest integer that fulfills n^L = 136.96 A, where 
L is the length of the ribbon in the periodic direction. 
This corresponds to a 32-fc-point scheme for 4.28 A, the 
dimension of the smallest unit cell of an armchair edge. 
We minimize the energy with respect to ionic position 
in the GNR plane; the relaxation is stopped when the 
force on each atom is less than 0.01 eV/A. We first de- 
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tcrminc the honeycomb lattice constant ao of the infinite 
graphene sheet, which is a = 2.471 A. This corresponds 
to a C-C bond length of ao/y/3 = 1.427 A. This agrees 
very well with the value of ao = 2.468 A found in another 
density-functional calculation^. 
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FIG. 1: A graphene edge orientation is defined by the angle 
a between edge tangent vector T e = nai + ma.2 and the 
armchair vector a ac = ai + a.2- The two highest symmetry 
edges are shown: the a = 0° (1,1) armchair edge and the 
a = 30° (1,0) zigzag edge. Also show is the more general (4, 1) 

graphene edge, corresponding to a = cos -1 ^^7f^ ~ 19.1°. 

A general graphene edge orientation can be described 
by a convention that is commonly used to describe the 
chirality of a single-wall carbon nanotube^. We use a 
to describe the angle between the tangent to the edge 
T e = nai + ma.2 and a second vector parallel to the 
armchair, as shown in Fig. [TJ This angle is given by 

a = cos -1 \J 4[{ n+m)™¥ ^n\- Naturally, the armchair and 
zigzag edges correspond to a = 0° and 30°, respectively. 
Because of the symmetry of the honeycomb lattice, we 
only need to consider angles within the range of 0° < 
a < 30°. In the present study, we study a discrete set 
of edge orientations in this a-range: 0°(1, 1), 6.59°(3, 2), 
10.89°(2,1), 16.10°(3,1), 19.11°(4, 1), 21.05°(5,1), and 
30°(1, 0). Figure^} shows a (5, 1) graphene nanoribbon, 
corresponding to a = 21.05°; each repeat period of this 
edge structure contains 5 zigzag units (Fig. [2p) and a 
kink that has the armchair structure (Fig. 

In the absence of hydrogen, the edge energy (per unit 
length) E cdgc is 

E cdgc = r-[BGNR(a, w) - n c Ei], (1) 

where £ , gnr(o^,u;) is the energy of an optimized GNR 
with edge orientation a and width w that contains nc 
carbon atoms, and E\ is the energy per carbon atom in 
the infinite, flat graphene sheet with the lattice parame- 
ter that minimizes the energy. The factor of two in Eq. [T] 
accounts for the fact that a GNR has two edges. In this 
work, we also consider the termination of all of the dan- 
gling bonds on the unreconstructed edges by hydrogen 
atoms. In this case, the edge energy is 




FIG. 2: (Color online) The edge structure of unreconstructed 
graphene sheets in the (a) armchair a = 0° (1,1), (b) a — 
21.05° (5,1), and (c) zigzag a = 30° (1,0) orientations. 



where £ , gnr+h(«,w) is the energy of the graphene 
nanoribbon where the dangling bonds associated with 
carbon atoms on the edge are saturated with single H 
atoms and njj is the total number of such H atoms. _Eh 2 
is the energy of the hydrogen molecule and the factor 
of 2 in the last term accounts for the fact that each H2 
molecule only contributes one H atom to each initially 
unsaturated bond. Note that Eq. [5] reduces to Eq. Q] in 
the absence of hydrogen tih = 0. 

By analogy with the surface stress^, we define the 
edge stress r c as the work required to stretch a graphene 
edge by an infinitesimal strain e (parallel to the edge), 
AE'odgo/e, or in differential form by 
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(Note: two definitions of surface stress t s have been 
widely employed in the literature: (a) t s = dj/de and 
(b) t s = 7 + dj /de, where 7 is the surface energy. We 
apply the approach used for surface stress definition (a) 
here, i.e., r e = dE e d ge /de. This is consistent with the 
method employed by Jur>i£, but not by Huang et alr^ in 
their studies of graphene edges. This definition yields the 
edge stress that produces elastic distortions of graphene 
nanoribbon edge a 15 ' 16 .) The derivative in Eq. [3] is deter- 
mined using a method based upon the bulk stress cal- 
culation within the Siesta code22 that assumes an affinc 
transformation of the atomic coordinates according to 
the applied strain. Throughout this work, we adopt the 
convention that a negative (positive) edge stress corre- 
sponds to compression (tension). 



III. EDGE ENERGY 

The first question we address, is how wide must a 
GNR be in order that the interactions between the op- 
posite, parallel edges arc negligible. Figure [3] shows 
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the edge energy versus GNR width, without and with 
H-termination. These figures clearly demonstrate that 
the edge energy converges to a nearly width-independent 
value for w > 30 A for all GNR orientations, both with- 
out and with H-termination. We note that earlier first 
principles have both included^ 7 - and exclude d 16 ' 26 spin 
polarization effects. We have not included spin polariza- 
tion effects here as they are expected to be more impor- 
tant for a proper determination of the band structure of 
zigzag edge graphene nanoribbons than for the overall 
edge thermodynamics. 
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FIG. 3: (Color online) The edge energy as a function of GNR 
width of the unreconstructed graphene nanoribbon (a) with- 
out hydrogen termination, and (b) with hydrogen termination 
for several GNR edge orientations. 

Unlike for all of the other edges, the a = 0° (1,1) 
armchair edge energy exhibits a decaying oscillatory be- 
havior with respect to GNR width w, as seen in Fig. [3] 
This is consistent with earlier observations for the spe- 
cial case of an armchair GNR without H-terminatio n 27 ' 28 . 
Son et alA demonstrated that the energy bandgap of 
armchair GNRs can be grouped into 3 families accord- 
ing to the GNR width, w. These families correspond to 
N a = 3p, 3p + 1 and Zp + 2, where p is an integer and 
N a is the number of carbon rows parallel to the armchair 
GNR axis. Figures 2] a and b show the edge energy ver- 
sus GNR width for the unterminated and H-terminated 
cases, where the data has been separated according to 
these families. With this division, the armchair edge en- 
ergy versus w largely converges monotonically. 

The converged (with respect to w) edge energies shown 
in Fig. [3] are summarized in Fig. [S^t, as a function of 
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FIG. 4: (Color online) The a = 0° (1,1) armchair edge as 
a function of ribbon width w grouped according to the di- 
vision of the data into three families^ (a) without hydrogen 
termination and (b) with hydrogen termination. 
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FIG. 5: (Color online) (a) The graphene edge energy as a 
function of the edge orientation a for the largest width GNRs 
in Fig. |3j The insets show the same data as in the main fig- 
ure but over a narrower edge energy range. The figure shows 
data for the unreconstructed edges both without and with 
H-termination. Also shown are data points for the recon- 
struction of the a = 30° (1,0) zigzag edges [ZZ(57)] without 
and with H-termination. (b) The edge stress as a function of 
edge orientation a for the same configurations as in (a). The 
hydrogenated zigzag(57) has too large of a tensile edge stress, 
2.570 eV/A to show on this plot. 
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edge type, a. For the unreconstructed and unterminatcd 
edges, this data shows that the edge energy increases 
nearly linearly with a; the armchair (zigzag) edge has the 
lowest (highest) edge energy, 1.202 eV/A (1.543 eV/A). 
The energy ordering (armchair versus zigzag) can be un- 
derstood by noting that the armchair edge reduces its 
energy by forming triple bonds in the armrest regions, as 
suggested by the fact that this bond length 1.258 A is 
much smaller that the bulk graphene C-C bond length of 
1.427 A. In contrast, the zigzag edge does not have the 
opportunity to form triple bonds; the outer C-C bonds 
relax slightly to yield a bond length of 1.400 A. The edge 
energy difference between the armchair and zigzag edges 
is 0.341 eV/A without reconstruction or H-termination, 
in agreement with previous DFT results^ (0.33 eV/A). 
The near linearity of the edge energy with respect to 
a can be understood by noting that all of the unrecon- 
structed edges are made up of a combination of armchair 
and zigzag units in proportion to the edge angle; i.e., 
steps on near armchair oriented edges are zigzag units 
and steps on near zigzag oriented edges are armchair 
units (see the (5,1) edge in Fig. 



IV. EDGE STRESS 

The graphene edge stresses, with and without hydro- 
gen termination, are presented in Fig. 0). Without hy- 
drogen termination, the edge stresses are compressive 
(negative) for all a. For the armchair [a = 0°, (1,1)] 
edge t c = —2.451 cV/A and for the zigzag [a = 30°, 
(1,0)] edge T e = -2.551 cV/A. These values are of the 
same order of magnitude as those obtained by Juni£ (i.e., 
armchair edge - r e = —2.640 eV/A and zigzag edge - 
t c = —2.248 eV/A). Calculations of the edge stresses us- 
ing an adaptive intermolecular reactive empirical bond 
order (AIREBO) potential, significantly underestimate 
the magnitude of the armchair and zigzag edge stresses 
(i.e., they find that the edge stresses of the armchair and 
zigzag edges are r e = —1.05 eV/A and r c = —2.05 eV/A, 
respectively). While earlier studies only considered the 
edge stresses of the two highest symmetry edges, the re- 
sults of Fig. 0) show the edge energy for many intervening 
orientations. Unlike for the edge energies, the edge stress 
is not a monotonic function of the edge orientation angle. 
Rather, the magnitude of the edge stress is a minimum at 
a ~ 15°. Reconstruction of the zigzag edge to from the 
zigzag(57) structure considerably raises the zigzag edge 
stress to near zero (slightly tensile). 

The main conclusion from this part of the study is 
that in the absence of hydrogen, all edge orientations ex- 
hibit compressive edge stresses (except for the zigzag edge 
upon reconstruction). We can understand the sign of the 
edge stress by consideration of edge relaxation (i.e., like 
surface relaxation) . Finnis and Heine^ presented an ar- 
gument in which surface atoms relax inward (toward the 
bulk) in order to smooth the electron density near the 
surface. This lowers the kinetic energy of the electrons, 



while having a relatively weak effect on the pseudopo- 
tcntial contribution to the total energy. The same thing 
happens at edges in the present two dimensional mate- 
rial. When the edge atoms move inward, they create a 
force on the atoms below the edge that is parallel to the 
edge and tends to push these atoms apart. However, the 
atom spacing, parallel to the edge, is fixed by the con- 
straint of the atoms in the bulk of the graphene sheet. 
It is these forces that create the compressive edge stress. 
When hydrogen adsorbs on the surface (see below), the 
inward relaxation effect is strongly reduced (due to near 
perfect saturation of the dangling bonds), thereby greatly 
reducing the magnitude of the compressive edge stresses. 

V. HYDROGEN ADSORPTION 

When hydrogen terminates all of the dangling bonds 
of the ribbon edges, the unreconstructed edge energies 
drop by 1 — 2 orders of magnitude, as seen in Figs. [SJa 
and [5^i. The edge energy in the presence of hydro- 
gen is maximum for the zigzag [a — 30°, (1,0)] edge 
(-Eedge = 0.103 eV/A) and minimum for the armchair 
[a = 0, (1,1)] edge (£ cdgc = 0.019 eV/A). While the 
magnitude of the edge energies drop dramatically upon 
hydrogen saturation, the ratio of the zigzag to armchair 
energies shows a much smaller change. 

Saturating the GNR edges with hydrogen significantly 
decreases the magnitude of the edge stresses (see Fig. [5b). 
Interestingly, saturating the zigzag GNR edge with hy- 
drogen changes the edge stress from highly compressive, 
-2.551 eV/A, to nearly zero, -0.006 eV/A. This is con- 
sistent with the results of an earlier study^ that show a 
significant drop in the magnitude of the edge stress upon 
saturating the edges with hydrogen. 

In order to understand the effect of hydrogen on 
the GNR edges, we focus first on the hydrogen ad- 
sorption energy £7 a( j s on the graphene edge: £7 a d s = 
(E-Gnr + 'Jh-Eh — £gnr+h) /ih- If we denote the num- 
ber of H atoms adsorbed per unit length of the edge 
as ?/h = nn/2L 7 and the binding energy of the hydro- 
gen molecule (relative to the isolated hydrogen atoms) 
as Eh = 2En — Eh 2 , we can rewrite i? c d go +H (sec Eq. |5J) 
as 

Eh 

-^edge+H = ^edge — ^H^ads — ). (4) 

This change in edge energy (per unit length) upon satura- 
tion with hydrogen is simply the product of the number 
of H atoms adsorbed (per unit length) and the change 
in energy per hydrogen atom in going from H 2 in the 
gas to H on the edge. For the unreconstructed edges of 
graphene, 77H has a very simple form: 

cos a /4 



where qq is the honeycomb lattice constant. ?/h is 15.5% 
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larger for the armchair (a = 0°) than for the zigzag edge 
(a = 30°). 
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FIG. 6: (Color online) The hydrogen adsorption energy as a 
function of angle a. The binding energy of a hydrogen atom 
in H2, Eb/2, is shown for reference. 

The main effect of adsorbing hydrogen on the edges is 
to decrease the edge energy by r](E a( i s — Figure^ 
shows that the energy of the zigzag edge decreases upon 
hydrogen adsorption more so than does the energy of the 
armchair edge. This is contrary to what we should ex- 
pect based upon the number of adsorbed hydrogen atoms 
per unit length. 77H (which decreases with increasing a as 
cos a). However, the hydrogen adsorption energy E a ds in- 
creases with increasing a (as seen in Fig.[f5]) more quickly 
than r?H decreases. This effect is associated with the 
fact that, before hydrogen termination, carbon atoms 
on the armchair edge have triple bonds to neighboring 
edge carbon atoms (and no dangling bonds), while the 
zigzag edge has no triple bonds (but docs have dangling 
bonds) 2 -. This suggests that hydrogen should bind more 
strongly to the zigzag edge atoms than to the armchair 
edge atoms. This implies (-Egnr — -E'gnr+h)/"-h will be 
larger for the zigzag edge than for the armchair edge. Be- 
cause the only term in -E ac i s that depends on edge type is 
(-Egnr — -E/gnr+h) /^h, this implies higher H adsorption 
energy to the zigzag edge than to the armchair edge. 

The hydrogen adsorption energy on the zigzag(57) 
edge, 5.22 eV/atom, is significantly smaller than the ad- 
sorption energy on the unreconstructed zigzag edge (see 
Fig. [5]). This can be attributed to the fact that triple 
bonds formed in the heptagonal units on the zigzag(57) 
edge and to the lack of triple bonds on the unrecon- 
structed zigzag edge. The difference in adsorption energy 
between the two zigzag edges explains why the hydro- 
genated zigzag(57) edge energy is (0.249 cV/A) higher 
than the hydrogenated unreconstructed zigzag edge (see 
Fig. 00 • It is important to focus upon H adsorption en- 
ergy -E a ds for the zigzag and hydrogenated zigzag edges 
(presented above) relative to the binding energy of hydro- 
gen in molecular hydrogen E^/2 (on a per atom basis). 
This comparison shows that while reconstruction reduces 
the tendency for H-passivation of the zigzag edge, H will 



still bind to this edge. (This conclusion differs from the 
interpretation of the discussion in Ref. l26l .) 

VI. EQUILIBRIUM GRAPHENE SHEET SHAPE 

In this section, we investigate the equilibrium shape 
of a flat graphene sheet, following a procedure akin to 
that used for determining the equilibrium shape of a 3- 
dimensional crystal via the Wulff construction (e.g., see 
Ref. 01). The first step is to construct a polar plot of 
the edge energy £ e dge as a function of the edge orien- 
tation angle a. This is shown in Fig. [7Jt for the case 
of pure, unreconstructed graphene edges using the data 
from Fig. [5^. Next, we construct a ray from the origin 
at edge orientation angle a on this polar plot. Then, we 
draw a line perpendicular to this ray at the point of its 
intersection with E e d ge - This process is repeated for ev- 
ery point (many points) around the polar E e dgc versus 
a plot. The inner envelope of these lines represents the 
equilibrium graphene flake shape. This construction and 
the resultant equilibrium flat graphene flake shape (i.e., 
that which minimizes the edge energy at fixed flake area) 
is shown in Fig. [7^, for the case of graphene with no edge 
reconstruction and without hydrogen termination. The 
equilibrium shape is a hexagon with flat (1,1) armchair 
edges. 
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FIG. 7: (Color online) Polar plots of the graphene edge energy 
_Eed gG as a function of a (a) for edges without hydrogen ter- 
mination or reconstruction, (b) with reconstruction but with- 
out hydrogen termination based upon a linear interpolation of 
the equilibrium armchair and zigzag(57) edge energies and (c) 
with hydrogen termination (no reconstruction occurs). The 
edge energy scales in (a) and (b) are identical, but are mul- 
tiplied by a factor of 16.7 in (c) for clarity. The red circles 
indicate data from our density functional calculations, the 
black lines are an interpolation of the data and the blue lines 
are drawn as perpendiculars to the radius -E e d g o vectors. The 
equilibrium, flat graphene flake shape is the inner envelope of 
the blue lines. 

The data presented in Fig. [5^ suggests that, in the ab- 
sence of hydrogen, the zigzag edge reconstructs and that 
the reconstructed surface has lower energy than even the 
armchair edge (which does not reconstruct). Since we 
do not have data for the reconstruction of all edges with 
0° < a < 30°, we make the simplifying assumption that 
the energies of these edges can be represented by a lin- 
ear interpolation between the energies of the unrecon- 
structed a = 0° armchair edge and that of the recon- 
structed a = 30° zigzag(57) edge (as holds for the un- 
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reconstructed edges - see Fig. [5^,). The resultant polar 
-^cdgc versus a plot is shown in Fig. [TJd. Unlike the po- 
lar -Eedge versus a plot for the unreconstructed case, this 
one is relatively smooth and nearly circular, yet a six-fold 
symmetry is obvious. Performing the Wulff construction 
on this interpolated E e dg e versus a plot produces an equi- 
librium graphcnc flake shape that has six-fold symmetry, 
but no flat surfaces. The absence of flat surfaces or facets 
can be traced to the fact that the polar -E e d g e versus a 
plot here shows no sharp cusps (discontinuity of slope). 
More interesting, however, is that when we allow edge 
reconstruction to occur (in the absence of hydrogen), the 
resultant six-fold equilibrium shape is rotated relative to 
that without edge reconstruction. That is, the equilib- 
rium graphene flake shape is dominated by surfaces at or 
near the zigzag edge orientation. 

When hydrogen is introduced into the system, such 
that all of the dangling bonds on the edges are saturated, 
the edge energy is greatly reduced, as discussed above. 
The hydrogenation of the edges removes any reconstruc- 
tion and the (1,1) armchair 0° (1,1) edge has the lowest 
edge energy and the zigzag 30° (1,0) edge has the high- 
est. The resultant polar plot of Sedge versus a is shown 
in Fig.[7k. This plot shows strong cusps at the armchair 
orientations. Performing the Wulff construction, we find 
that the equilibrium graphene flake shape is hexagonal 
with (1,1) armchair edges, as in the unreconstructed, un- 
hydrogenated graphene edge case. 

The above analysis focused on the equilibrium shape 
of a flat, finite graphene sheet (a two-dimensional intrin- 
sic manifold). However, a flat graphene flake is unstable 
against thermal rippling at finite temperature^ .. Com- 
pressive edge stresses can also lead to elastic distortions 
of/near the graphene flake/ribbon edges, as discussed in 
the literatur o 13 ' 15 i 31 . These distortions may be described 
as the buckling of the edges, associated with the com- 
pressive edge stress but constrained by the elasticity of 
the interior of the graphcnc flake. The amplitude of the 
buckling A has been described^ as 
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FIG. 8: (Color online) The edge buckling amplitude for the 
pure and hydrogenated, unreconstructed edges, as per Eq. [6] 
using the data from Fig. [5b. 

magnitude of the buckling amplitude reported in Ref . fl5L 
The buckling amplitude is in the range of 2.70 to 3.11 A 
over the entire range of a for the unreconstructed edges 
without hydrogen. Interestingly, if the zigzag edge recon- 
structs to the (57) structure, the buckling amplitude, as 
per Eq. ([5]), is imaginary (the edge stress is positive); 
implying that the reconstructed zigzag edge does not 
buckle. 

When the edges are saturated with hydrogen, the edge 
stress drops dramatically (by a factor of ~ 7.54 to 425 
for the armchair and zigzag edges, respectively) and the 
buckling amplitude drops to 1.11 and 0.14 A for the arm- 
chair and zigzag edges (see Fig. This represents a 
drop in edge buckling amplitude by a factor of ~ 2.3 to 
22. In short, edge hydrogenation strongly reduces the 
edge stress driven graphene edge buckling to the point 
that it is, for most purposes, negligible. 



VII. CONCLUSIONS 



V (tjV 20 + 14y7/18)M / 

where A is the wavelength of the buckling, M = 1&/(1 — 
v 2 ), Yb and v are the Young's modulus and the Poisson 
ratio of the graphene sheet, and where we have neglected 
the relatively minor effect of the elastic stiffness of the 
edge. Assuming a typical wavelength of A = 10 rnn and 
M = 2000 eV/nm 2 for graphene^, we can use the data 
presented above to determine the buckling amplitude A 
as a function of edge orientation a and examine the effect 
of hydrogenation of the edge. 

Figure [5] shows the buckling amplitude from Eq. © 
using the edge stress data from Fig. [5b as a function of 
a, with and without hydrogen on the edges. The edge 
buckling amplitude found here is 3.05 and 3.11 A for the 
unreconstructed armchair and zigzag edges without hy- 
drogen, respectively, which is roughly consistent with the 



We have employed density functional theory-based cal- 
culations of graphene nanoribbons to determine the edge 
orientation dependence of the edge energy and edge stress 
(akin to surface stress in three-dimensional materials). 
In the absence of reconstruction, the edge energy is a 
nearly linear function of edge orientation angle; increas- 
ing from the armchair orientation to the zigzag orienta- 
tion. Reconstruction of the zigzag edge, however, lowers 
its energy to less than that of the armchair edge. The 
edge stress for all edge orientations is compressive. The 
magnitude of the edge stress is a minimum at an orien- 
tation nearly midway between the armchair and zigzag 
orientations. Reconstruction of the zigzag edge reduces 
the edge stress to near zero (very slightly tensile). Hy- 
drogen adsorption to graphene edges is favorable for all 
edge orientations, with a larger adsorption energy to the 
unreconstructed zigzag edge than to the armchair. Hy- 
drogen adsorption dramatically lowers all edge energies 
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and all edge stresses. It also removes the reconstruction 
of the zigzag edge. 

The thermodynamic properties of graphenc edges play 
a key role in determining the equilibrium shape of a finite 
graphene sheet or flake. Using the new edge energy data, 
we determined that the equilibrium shape of a graphene 
sheet (with unreconstructed edges) is hexagonal with 
straight armchair edges. This result is valid for both the 
hydrogen- free case with unreconstructed edges and in the 
presence of hydrogen. However, if the zigzag edges recon- 
struct (as thermodynamically favored), graphene flakes 
will have a six-fold symmetry, but with rounded (rather 
than faceted) edges. Interestingly, the six-fold symmetry 
is rotated with respect to the hydrogenated case and will 
be dominated by near zigzag edges. 

The compressive edge stresses will lead to edge buck- 
ling (out-of-thc-planc of the graphene sheet) for all edge 



orientations, in the absence of hydrogen. The edge buck- 
ling amplitude is expected to be approximately 3 A, but 
this value depends on buckling wavelength. Exposing the 
graphene flake to hydrogen will dramatically decrease the 
buckling amplitude to the point that it may be too small 
to be of any significance. 
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